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O RELATIVISTIC LIDOV-KOZAI RESONANCE IN BINARIES 

^ Cezary Migaszewski ' and Krzysztof Gozdziewski ' 

O 

Abstract. We consider the secular dynamics of a binary and a planet in terms 
7-H of non-restricted, hierarchical three-body problem, including the general rela- 

tivity corrections to the Newtonian gravity. We determine regions in the param- 
eter space where the relativistic corrections may be important for the long-term 
'^j ' dynamics. We try to constrain the inclinations of putative Jovian planets in 

recently announced binary systems of HD 4113 and HD 156846. 

Q-{ 1 Introduction 

6 

;h In the recent sample of detected extrasolar planetary systems, some planets exhibit large 

eccentricities. It may be explained by the Lidov-Kozai resonance (LKR) acting in binary 
^ ^ stellar systems (e.g.. |Innanen etal. 199 7 Takeda and Rasio 2005 Verrier and Evans 2008]). 

If the inner planetary orbit is inclined to the orbital plane of the binary, the exchange of 

I the angular momentum between orbits may force large amplitude eccentricity oscillations 

of the planetary orbit, and simultaneously its argument of pericenter coi librates around 
±7t/2. However, the LKR may be suppressed by the general relativity (GR) correction to 
the Newtonian gravity (NG) through changing frequencies of pericenters. Here, we focus 
on the non-restricted problem and relatively compact systems, and the dynamical effects 

• of including the GR interactions in the model of motion. 

o 

^ 2 The secular dynamics of the hierarchical triple system 



We consider the hierarchical triple system. The Hamiltonian written with respect to 
canonical Poincare variables (e.g.. lLaskar and Robutel 19951 1. y-f = i^epi + i^ert, where 



X 

f~^i V2p/ n J V ^ "^0 ) 

describes perturbed Keplerian motions of the inner binary (the central mass otq and m\), 
and the outer binary (mo and more distant point-mass m-i), /j* — (otq + m,), where k is 
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the Gauss gravitational constant, P,- = (l/m,- + are the reduced masses, ri.2, are 

the radius vectors of mi 2 relative to mo, pi,2 stand for their conjugate momenta relative 
to the barycenter, and A — ||ri — r2||. stands for GR correction to the Newtonian 



potential of wiq and m\ (see, e.g., Richardson and Kelly 1988) l. We assume that the ratio 
of semi-major axes a = a\/a2 < 0.2, and i^^ei-t ^ ^H^epi- It means that both mi. 2 are small 
(planetary regime) or one of mi. 2 ~ mo is relatively large, and one of these bodies is 
enough distant from mo (binary regime). 

We expand i?^G with respect to a and the Hamiltonian is averaged out with respect 
to the mean longitudes ( [Migaszewski and Gozdziewski 2008^ , that leads to the secular 
term {^g) + {H^v), where 



k m\m2 

fl2 



1 + J 1 - e2 ^ X'5?i(ei ,e2, (Oi , CO2,/) 



1=2 



X = a/{\-el), (2.2) 



/ stands for the mutual inclination, coi 2 are the pericenter arguments, and perturbing 
terms ^ are derived in (Migaszewski & Gozdziewski, in preparation). The averaged GR 
term is (^^^r) = — 3pijUjC^^aj"^(l — ej)^'/^, where c is the velocity of light. The expan- 



sion in Eq. 2.2 generalizes the octupole theory (e.g., IFord et al. 20001 ) and the coplanar 
model (Migaszewski and Gozdziewski 2008a I. After the Jacobi's elimination of nodes 
(Ail = ±Jt), we eUminate one degree of freedom thanks to the integral of the total angular 
momentum, C. Then i?Z|ec = ^ec(Gi , G2, COi , CO2) parameterized by C = |C| (or the Angu- 
lar Momentum Deficit AMD = Li +L2 — C, where L1.2, Gi 2 are the Delaunay actions) is 
reduced to two degrees of freedom. For a = 0.1, the relative errors of approximated 
by the 10-th order expansion do not exceed 10^^ in the relative magnitude (see Fig. [I}. 

To study O-Qec, we apply the representative plane of initial conditions, E, introduced 
in jMichtchenko and Malhotra"^00 4, M ichtchenko et al. 200 6 ) which crosses all phase- 
space trajectories. Due to symmetries of i^^ec with respect to the apsidal and nodal lines: 



3a)i 



3(02 



:0, (coi,co2)e{(0,0),(0,±7t),(±7t/2,±7t/2),(±7t/2,T7t/2)}, (2.3) 



and these conditions define the E-plane, E — {ei cosA!Ii,e2Cos2a)i}, ACJ = CJi — 052, 
(ei,e2) G [0, 1), see the left-hand panel of Fig. |2] for an illustration. Restricting (0)1,0)2) 
to the above set, we also define E5 = {ei sin 0)1,^2 sinO)2}, Ec = {ei coso)i,e2COSO)2} 
revealing levels of y-Qsc without discontinuities dLibert and Henrard 20071 ). 



3 Stationary solutions and the Lidov-Kozai resonance 

The equilibria of i?^ec provide much information on the structure of the phase space. In 
the E-planes, these equilibria appear as quasi-elliptic or quasi-hyperbolic (saddle) points 
of the levels of J-Qec, according with the equations of motion: 

^—==0, ——=0, or ^ — = 0, ^ — =0. (3.1) 
oOl 0C72 ae\ ae2 

The stability and bifurcations of equilibria in the full and in the restricted three-body prob- 



lem were studied in many works (see, e.g., iKozai 1962.,Krasinsky 1972 Krasinsky 1974 
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Fig. 1. A test of the relative accuracy of the 10-th order expansion of ^^c, and levels of i?4ec in the 
Z^-plane (see the text for details) for different values of AMD compared with the semi-analytical 
(exact) averaging (see lMichtchenko and Malhotra 20041 or |Migaszewski and Gozdziewski 2008b| >. 
Differences between the theories are expressed in terms of the relative log-scale. 



Lidovand Ziglin 1974 Fejoz 2002', ' Michtchenko et al. 20061 ILibert and Henrard 20071 
Migaszewski and Gozdziewski 2008b i regarding the NG model. Here, we investigate 
more closely the equilibrium at the origin (ei = e2 = 0), which is well known since 
Poincare, in the presence of the GR interactions. According with the terminology of 
|Krasinsky 1974, that is the trivial space solution of the 3rd kind (e — Q,I^ 0), see Fig.[T^. 
The zero-eccentricity equilibrium (ZEE) is related to the maximum of y-Q^^ and is Lya- 
punov stable. For a given value of C, the mutual inclination of circular orbits, /q, is also a 
maximal mutual inclination if I\ 2 < 7i/2. Moreover, for some smaller C (larger AMD), the 
origin may change its stability due to bifurcations illustrated in the E^-plane (Figs.[T|3,c). 
For instance. Fig. [TJ) illustrates a saddle accompanied by two elliptic points. Close to 
these points, the phase-space trajectories exhibit librations of C0i,2 around ±7t/2. This 
structure (see also Fig. [2]i is associated with the LKR; the elliptic points may be called 
nontrivial, negative solutions of the 3rd kind, (e ^ 0,/ ^ 0), as in (Krasinsky 1974). They 
appear when C < Cent, or, equivalently, when io > icm for initially circular orbits. For 
more details see e.g., (ILibert and Henrard 20071 Migaszewski and Gozdziewski 2008b i. 

Here, we restrict our calculations to i„it < 7i/2 (the case of direct orbits), hence we do 
not follow the second bifurcation (Fig.[T|;) appearing for i^m ^ ll/2,e\ ^ 1. We compute 
i'crit causing the stability change of ZEE in the NG model for mass ratio fj = mi /m2 G 
[10"^ 10^], and a e [10 ~^0.2] (see the right-hand panel of Fig. |2}. Two kinds of LK 
bifurcation may appear (Krasinsky 1972 1: at /crit ^ 40° (the inner LKR, amplifying ei) 
and for i^rit ^ 64° (we call it the outer LKR; e.g, in a case of a circumbinary planet). 
Moreover, in the planetary regime of mi 2, iait depends only on a and /j, and not on 
individual semi-major axes nor masses. 



4 Effects of the General Relativity correction 



After introducing the i?^}R correction to i?4ec, the structure of the phase space changes 
qualitatively (Fig.[3]i. We choose the same io (a function of constant Li 2 and C) for fixed 
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Fig. 2. The left-hand panel: the representative plane of initial conditions, Z. The right-hand panel: 
the critical inclination Jcrit in the (/j, a)-plane, the NG model. See the text for more details. 
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Fig. 3. The Z.s-plane for mo = Im©, a\ = 0.5au, a = 0.01, /j = 0.01, and iq = 60°. The left-hand 
panel is for the NG model, next panels are for the GR model, and mj = {l,0.1}mj, respectively. 



a = 0.01 and /j = 0.01. For the NG model, a clear LKR structure appears (Fig. [3^). 
However, in terms of the GR model, the saddle structure may shrink (Fig.jS]?), and finally, 
for small enough masses, it disappears (Fig.|3j;). That effect may be characterized globally 
through the critical inclination idt = /critl/JjOc) for varying m\,a\ (Fig.|4]). The structure 
of the (/J, a)-plane in terms of the GR model is very different from the NG case (Fig. |2|l. 
We may see three distinct regions related to the inner LKR (/crit ^ 40°), and to the outer 
LKR (i'crit ~ 64°), smoothly passing into a new region emerging in the bottom-left corner, 
which is colored in yellow, where /crit — * 7t/2, and the LKR may be totally suppressed. 

5 An application to the HD 4113 and HD 156846 planetary systems 

We apply the results to test a hypothesis that highly eccentric orbits of recently detected 
Jovian planets in HD 41 13 and HD 156846 planetary systems dTamuz et al. 2 008) can be 
explained by the LKH resonance (in the sense considered here, i.e., of the direct orbits) 
forced by more distant and unseen (likely massive) objects. Indeed, the radial velocity 
(RV) of HD 4113 exhibits annual trend RV; ~ 28 ms^^ that implies the minimal mass 
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Fig. 4. The critical inclination I'crit in the (jU, a)-plane, in terms of the relativistic model. Orbital 
parameters of the innermost body are labeled in the respective plots, mo = Img. 



of putative distant companion ~ lOmj and ~ lOau. A simulation of /crit is illustrated 
in Fig. [5^. Here, we assume that the orbit of HD 411 3b is edge-on. The red lines in the 
{fj, a)-plane mark raw limits of the orbital parameters of a putative distant object. It must 
be also massive enough to induce the observed RV drift, hence the skew line is for the RV, 
estimated under an assumption that its orbit is circular. Moreover, to force maxeb ^ 0.9 
(to generate large enough "loop" in the E-plane, see Fig.[3|l, appropriately large /q < 7i/2, 
('0 > 'crit), is required. Figs.|5|5,c illustrate min/o(maxe), i.e., the minimal inclination /q 
for which ei may reach given maxe (here, maxe = 0.9). Figure[5]3 is for the NG model 
and Fig. |5j; is for the GR model, respectively. In both cases, min/o(0.9) ^ 70° and the 
putative body may be responsible for the detected large eccentricity of HD 4 11 3b. 

The Jovian planet HD 156846b belongs to a wide binary with a > 250au and Wb ~ 
O.56m0. We found that the putative system is located in such a {/j, a)-region, in which the 
LKR can be suppressed at all because iait ^ 1^/2 (Fig.jsjl). Moreover, min/o(0.85) ^ 66° 
for the NG model (Fig.|5^), while minio(0.85) > %/2 for the GR model (Fig.|5f), also the 
structure of the a) -plane is qualitatively different in these two cases. Hence, in the GR 
model, the eccentricity cannot be explained by the LKR (in the sense considered here). 



6 Conclusions 



Recently (Migaszewski and Gozdziewski 2009 1, we found that apparently subtle GR cor- 
rection to the Newtonian model of coplanar planetary system may lead to significant, 
qualitative changes of the secular dynamics. In the present work, we try to extend such 
a quasi-global study to non-coplanar model, applying the averaging and the concept of 
representative plane of initial conditions. The results indicate that the 3D dynamics are 
also very different in the both 3D models. We continue the work on this problem. 
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Fig. 5. The critical inclination i^yit (the left-hand cohimn), min;o(maxe) for the NG-model (the 
middle column), and min!o(maxe) for the GR-model (the right-hand column). The top row is for 
the HD 4133 system (maxe = 0.9), the bottom row is for the HD 156846 system (maxe = 0.85). 
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